Charge-transfer polaron induced negative differential resistance and giant 
magnetoresistance in organic spintronics: A Su-Schrieffer-Heeger model study 
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Combining the Su-Schrieffer-Heeger model and the non-equilibrium Green's function formalism, 
we investigate the negative differential resistance effect in organic spintronics at low temperature 
and interprete it with a self-doping picture. A giant negative magnetoresistance exceeding 300% is 
theoretically predicted as the results of the negative differential resistance effects. 
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The discovery of negative differential resistance (NDR) 
in traditional semiconductor diodes and also organic 
semiconductor (OSE) nanostructures 0, El has opened 
a new chapter of device physics. Motivated by the po- 
tential applications of organic NDR, numerous experi- 
ments have been done with different OSE electronic de- 
vices 0, 0, 0) 01- Although some possible mechanisms 
have been suggested, the organic NDR remains a theo- 
retical challenge, beyond the simple picture of interband 
tunneling or resonant tunneling in heterostructures. The 
organic NDR concerns not only the excess c harg e transfer 
through lead/OSE/lead structures 0, IE 0> E| > but also 
the strong electron-phonon (e-vh) coupling that induces 
polarons in OSE structures (lli . 

Historically, the discovery of magnetoresistance 
(GMR) effect in 1988 is considered as the beginingof 
a new technology called spintronics j^J [Tj, [l!| Hs LLd] , 
where it is not only electron charge but electron spin that 
carries information. In addition to inorganic semicon- 
ductors 0, El ) OSE materials, due to their controllable 
structure, strong e-ph coupling [Tl|, and large spin coher- 
ence |20|, offer another promising system to spintronics 
|2ll |22| |. In a recent experiment of spin-injection in a 
spin- valve structure LSMO / A\q^ /Co, a GMR as large as 
40% had been detected 22]. In this Letter, we report 
the possibility of a charge-transfer polaron-induced NDR 
mechanism via a numerical study on a representing tt- 
conjugated OSE model device that exhibits also a giant 
GMR. 

In the study of the electronic conductivity and opti- 
cal phenomena in tt conjugated OSEs, the Su-Schrieffer- 
Heeger (SSH) model 23] has also been shown a remark- 
able track of success. This model captures the essen- 
tial characteristic of a conjugated molecule, where the 
strong e-ph coupling leads it to the polaron (or solition) 
charged states and dimerized ground state. In addition 
to polyacetylenes, the SSH model has been applied to 
charged conjugated systems [24|, carbon nanotubes plj . 



and DNA molecules [2J, |2J. In this work, the SSH 
Hamiltonian for the OSE electrons coupled adiabatically 
with lattice displacements reads 
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Each atomic unit in OSE is represented by a single nor- 
malized site; c+ a {c nM ) denotes the creation (annihila- 
tion) operator of an electron at the n th site with spin 
a, while e a , t a , and t\ are the on-site energy, zero- 
displacement hopping integral, and nondegeneracy pa- 
rameter, respectively. The lattice distortion is treated 
classically, in terms of the bond distances {y n — it„+i — 
u n }, deviated from its energy- minimum values that are 
to be determined via the Hellman-Feynman theorem [cf. 
Eq. JfjJ] , with the spring constant K Q and the adiabatic e- 
ph coupling constant a . This is a static polaron model, 
which together with the effective noninteracting many- 
electron ansatz in Eq. QJ, are considered to be justifiable 
in the present study of stationary transport. The effects 
of polaron motion and electron-electron correlation (be- 
yond the Hartee-Fock approximation) will be the subjects 
of forthcoming work. 

Further, we choose a symmetric ferrimagnetic (FM) 
3c? transition metal as electrodes. The spin-dependent 
charge transport takes place between their 3d bands. 
Neglecting spin-flip during transport and adopting the 
two-current model 28], we describe the FM metal by 
one-dimensional single d-band tight-binding model with 
a spin splitting term pd| , 

H F = ^|e / d+ CT d„, . + t/(d+ CT d n+ i >0 . + H.c.)| 

n,cr 

-E J /(<At-<Ai). ( 2 ) 
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where a (d n _ a ) is the creation (annihilation) operator 
of an electron in the metal at the n th site with spin a; 
tf is the on-site energy of a metal atom, tf is the nearest 
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neighbor transfer integral, and J/ is the Stoner-like ex- 
change integral. The coupling between the OSE and FM 
electrodes is described by the spin-independent hopping 
integral, £ c l = £ c r = 0(tf + to), where (3 denotes the 
OSE-metal binding parameter. 

The non-equilibrium Green's function (NEGF) ap- 
proach based on the Keldysh formalism [2|| [23, 0HJ is 
used to calculate the quantum transport properties of 
organic spintronics. To do that, the spintronic device is 
divided into three distinct regions. One is the so-called 
central scattering region (S-region), with the Hamilto- 
nian H$ = i?L + Ho + Hr + H int , which consists of 
the OSE together with a small number of metal atoms 
attached to each of its ends. The other two are elec- 
trodes (L and R) that serve as charge reservoirs with 
the steady state electronic distribution of bulk metal at 
given temperature. Tracing out the reservior degrees of 
freedom leads to an effective S-region Green's function 
G{E) = [ES -H s - E L (£) - E R (E)] -1 - Here, S (set 
to be the unit matrix) is the overlap integral matrix be- 
tween basis wave functions, while E L / R is the self-energy 
matrix that accounts for the effects of reservior electrodes 
on the S-region [33] • It is possible to have the analyti- 
cal solution of E L / R for the one-dimensional FM metal 
transfer-coupling with the OSE system. In this work, we 
adopt an efficient numerical approach throu gh s olving 
eigenvalue equations to achieve the self-energy [3_2j . This 
approach can be easily extended to study the magneto- 
transport beyond the one-dimensional system. 

The current can now be evaluated as [30l | 

oo 

I= T l Tr ( r L Gr « Gt )[/( S ^L)-/(i?,/iR)]^. (3) 

— oo 

The trace term in the integrand is the transmission co- 
efficient function, in which Tl/r = ^(El/r — ^l/ r ) de- 
notes the broadening matrix; f(E, /^l/r) is the Fermi dis- 
tribution function at the lead chemical potential /^l/r- 
One can also evaluate the density of states (DOS) via 
D(E) = iTr[G(E) - G t ( J B)]/(2 7 r), and the reduced den- 
sity matrix as 

oo 

P=Y J2 [ GT a G^f(E,fi a )dE 7 (4) 

with a given number of carrier electrons N = Trp in the 
S-region. We separate the reduced density matrix into its 
equilibrium and bias-induced contributions, and evaluate 
them by contour integration and direct multi-grid Gaus- 
sian integration, respectively. Numerical implementation 
will be carried out in a real-function basis set represen- 
tation, and the resulting p will be real [3l|. 

We take the bias V not changing the electronic struc- 
tures of L and R reservoirs, but just shifting their poten- 
tials by V/2 and — V/2, respectively. In contrast, it does 
alter the Hamiltonian of OSE to become Ho — ecj>(x). 
The electric potential drop here, satisfying <fi(xi) = V/2 



and 4>(xn ) = —V/2, should be evaluated via Poisson's 
equation with the method of images j33| 

V 2 cp{x) = -Y j ^8{x n -x) 1 (5) 

with the charge density p n ^ n that depends also on (j>(x), 
the vacuum permittivity e n , and the n th site coordinate 
x n = (n— l)a, where a is the OSE lattice constant. 

Furthermore, the OSE lattice distortion, which is cor- 
related with electron wavefunctions due to the strong 
e-ph interactions, is assumed to have catched up with 
the charge variation without draggling. The resulting 
lattice distortion at the finite bias voltage can there- 
fore be evaluated via the Hellman-Feynman theorem: 
d[Tr(H s p)}/du n = 0; i.e. [cf. Eq. Q], 

2a (p„,„-i - Pn+i, n ) + K (y n -i - y n ) = 0, (6) 

with the index n including only the OSE sites, since there 
are no lattice distortion in the L and R sub-structures. 
From Eqs. (@J-@), one can see that both the charging 
effect from the electrodes and the external potential from 
the bias voltage are all included in the non-equilibrium 
density matrix of the coupled e-ph system which should 
be evaluated self-consistently. 

There are two distinct transport measurement config- 
urations, parallel (P) and antiparallel (AP), with respect 
to the relative magnetization orientation of two FM elec- 
trodes. Generally, in the two-current model, both the 
majority-majority and minority-minority (or majority- 
minority and minority-majority) transports are permit- 
ted in the P (or AP) configuration. But for the cobalt 
electrodes in the present work, the full-filled up-spin 
electrons cannot transport under in the P-configuration 
and the charge carriers are only the half-filled down- 
spin electrons. Contrary, the charge carriers in the AP- 
configuration under V > are only the up-spin elec- 
trons, driven from the full-filled majority-spin band of 
L-electrode to the half-filled minority-spin band of R- 
electrode of opposite magnitization orientation. 

We are now in the position to elucidate numerically 
the I-V characteratics, especially the NDR behavior of 
OSE systems. The SSH parameters for the OSE system 
[Eq. (JTJ] are e a = -4.3 eV, t Q = 2.5 eV, h = 0.04 eV, 
a Q = 5.0 eV/A, and K a = 21.0 eV/A 2 , which are the 
modified values from those of conducting polymers [34[ 
for the OSE being of larger band gap. The tight-binding 
parameters for the cobalt electrodes [Eq. are e/ = 
—7.0 eV, tf = 1.5 eV and J/ = 1.45 eV, determined by 
fitting the cobalt 3d band structure [35J. The OSE-metal 
binding parameter (3 will be specified later. 

To justify the above values, we calculated the band 
structures of the model OSE and FM cobalt electrodes 
individually (with /3 = 0); see Fig.^, for their result- 
ing DOS at T = 11 K. Our model OSE of 10 sites has 
the highest occupied and lowest unoccupied molecular or- 
bitals of -Ehomo = E F - 0.9 eV and Blumo = E F + 2.0 
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eV, where Ep is the Fermi level of the cobalt electrode. 
These results well reproduce the band gap diagram of the 
OSE spin- valve device used in a recent experiment |2^ . 
The model cobalt metal is also consistent with the real 
material |35j, which shows one full- filled majority-spin 
band, and one half-filled minority-spin band. 

Let us start with the equilibrium property of the 
model spintronic system at V = 0, at which the S- 
rcgion (Co/OSE/Co) assumes charge neutral. In our 
model, the S-region consists of a cobalt oligomer of 10- 
sites to each terminal of the 10-site OSE; thus each seg- 
ment in Co/OSE/Co involves 20 spin-orbitals. If the 
OSE-metal binding parameter /3 = 0, three segments of 
the S-region would be charge neutral individually, with 
A^l = =15 and Nq = 10 electrons, respectively. In 
reality /3 ^ and intra-regional charge transfer (ICT) is 
possible. Depicted in Fig.^}) is the calculated DOS of the 
organic spintronic device with the OSE-metal binding pa- 
rameter (3 = 0.5, where an ICT of about 1.14e from OSE 
to Co-segments has occured. The observed 'self-doping' 
phenomenon here is largely due to the strong e-ph inter- 
action, which leads to the formation of a preexisting hole 
polaron that stabilizes the S-region complex before ap- 
plying potential bias |36|. The preexisting polaron state 
is rather evident by examining the majority-spin band 
of Co/OSE/Co complex (solid-curve in Fig.QJ)), since its 
isolated metal counterpart (thick-solid curve in Fig.^,) 
is completely filled up to the Fermi level. 

We then calculated the I-V characteristic [Eq. ©] of 
the model Co/OSE/Co spintronics at T = 11 K and 
(3 = 0.5, in both the P and AP configurations of rela- 
tive magnitization orientation of FM electrodes. The re- 
sulting I-V curves in these two configurations (P: solid; 
AP: dash) are shown in Fig.[2jL, and the corresponding 
dl/dV ones are in Fig.[3). Included in Fig.|2b is also 
a thick-curve for the bias voltage dependence of magne- 
toresistance (MR), AR/R = (Rap — Rp)/Rap, measur- 
ing the relative difference of electric resistance with these 
two configurations. 

Consider first the NDR behavior (about — 26.3 kfi at 
its minimum) in the P-configuration, where the current, 
after an initial near-ohmic increase, drops quickly from 
/peak = 4.45 ^A at V poak = 0.35 V to I vaUey = 0.49 //A 
at V va \i C y = 0.4 V. To see what happens during the NDR 
region, we also examined other noncquilibrium proper- 
ties (at V > 0) . Shown in Fig.0 are the representing 
results of both the current peak (solid) and valley (dash) 
states: (a) the majority-spin (up-spin) DOS D(E); (b) 
the P-configuration (down-spin) transmission coefficient 
function T(E). Indicated in Fig .[3 are also the numbers 
of electrons in the OSE segment at the two corresponding 
voltages. By checking the charge distribution and the lat- 
tice distortion (not shown here), we found: (i) The pre- 
exising hole polaron remains localized around the OSE 
center when < V < V poa k; (ii) As V increases further, 
the excess electron charge migrates from the leads into 
the OSE segment; and (iii) At Kaiioy = 0.4 V, the pre- 
existing hole is completely annihilated, and the OSE is 



essentially in its dimerized ground charge-neutral state. 

The above observations suggest that the NDR be- 
haviour in the P-configuration, shown by the solid curves 
in Fig.|3 is due to the annihilation of the preexisting, 
'self-doping' polaron. As discussed earlier (cf. Fig.njj), 
the preexisting hole polaron level at V — deeply local- 
izes in the down-spin band (P-configuration conduction 
band) gap of the OSE, which leads to its relatively large 
DOS, and thus a low-resistance state according to the 
doping theory of conducting polymers jj^J ■ This accounts 
for the rapid increase of current when V < Vp ea k in the P- 
configuration (the solid curve in Fig.|2Ji). Figure [3] shows 
clearly that the increase of bias voltage from V pca k to 
Vvaiiey accompanies with the annihilation of the preexist- 
ing polaron, which accounts for the NDR observed in the 
P-configuration (solid curves) in Fig.|3 

In the AP configuration, the charge carriers at V > 
are no longer the down-spin electrons, but the up-spin 
ones, from the majority subband of L to the minority 
subband of R electrode. As the observed switch-on volt- 
age in this case (about 0.5 V) exceeds the aforementioned 
NDR region, the 'preexisting polaron' makes no direct 
contribution to the conductance in the AP configuration. 
The resulting I-V characteristic (dash-curve in Fig.|2i) 
can thus be all understood (including its switch-on volt- 
age) by examining the structures of the two involving 
subbands in Fig.[TJ) at V — 0. 

Finally, let us make comments on the voltage- 
dependent MR (the thick-curve in Fig.^a), especially 
the negative GMR of (-AR/R) max = 300% at V = 
1.1 V. Negative MR has been experimentally observed 
in the Co/SrTi0 3 /LSMO tunnel junction [33 and the 
LSMO/OSE/Co spin- valve device ^ . Traditionally, 
one analyses the observed MR, or other transport behav- 
ior as function of bias potential, via the involving DOS 
of conduction bands/subbands of uncorrelated ((3 = 0) 
FM/OSE/FM systems at V = 0, such as Fig.[TJi and 
Ref. l35l This traditional analysis does lead to our un- 
derstanding the qualitative MR-V behavior in each indi- 
vidual voltage range in Fig.03, where AR/R decreases 
from 100%, changes sign into the negative MR region, 
reaches (—AR/R) max at V = 1.1 V, and etc. However, 
the quantitive GMR values in Fig.|3), especially its max- 
imum value of 300%, cannot be accounted for via the 
simple band structure analysis with [3 = and/or V = 0. 
The calculated extraordinary GMR of 300% can only be 
accounted for via the 'self-doping' of preexisting polaron 
and its annihilation that affect distinctly differently on 
conductances depending on the (P or AP) relative mag- 
netization orientation of FM electrodes. 

In summary, we proposed to exploit 'self-doping' 
FM1/OSE/FM2 structures for distinct NDR and GMR 
materials. Numerical demonstrations were performed 
based on a realistic model system. The theoretical NDR 
was found as a result of transition from the low-resistant 
preexisting hole (cation) polaron to high-resistent dimer- 
ized charge-neutral ground state. This NDR mechanism 
is different from the two-step reduction picture, proposed 
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originally by Reed and Tour and co-workers [2( in ex- 
plaining their experiment on a redox-center-containing 
molecule, in which NDR results as transition from con- 
ducting anion to insulating dianion state. As it appears 
in the P- but not AP-configuration of relative magnetiza- 
tion orientation of FM electrodes, the NDR leads also to 
a large magnitude (300%) of negative GMR, much larger 
than the experimentally reported ones by far |22( ■ Many 
OSEs are easy-doping materials j34l |. and the required 



'preexisting polaron' state can be formed in cither self- 
doping or external doping manner. Thus, experimental 
realizations of GMR far exceeding 100% in organic spin- 
tronic devices should be feasible according to the present 
model study. 
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FIG. 1: DOS of the majority-spin (solid-curves) and 
minority-spin (dash-curves) bands of model systems: (a) The 
isolated OSE (thick-curves) and Co 3d band (thin-curves), 
with the OSE-metal binding parameter (3 = 0; (b) The binded 
Co/OSE/Co spintronic device with (3 = 0.5. Temperature 

r = iiK. 



FIG. 2: (a) The I-V characteristics of the Co/OSE/Co sys- 
tem in Fig.^i, measured with the P- (solid-curve) and AP- 
configuration (dash-curve) of relative magnetic orientation of 
electrodes, (b) The dl /dV, obtained numerically via (a), and 
—AR/R (thick-line) as functions of the bias voltage. 
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FIG. 3: The current-peak state (solid-lines) and the valley 
state (dash-lines), referring to the I-V charecteristic of the 
parallel configuration in Fig.|5^,, in terms of (a) the DOS of 
majority-spin band, and (b) the transmission coefficient of 
electrons (minority-spin). The total numbers of electrons in 
the OSE segment in the current-peak and valley states are 
specified. 



